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Abstract 

Planetary magnetic fields could impact the evolution of planetary atmospheres and have a role in the 
determination of the required conditions for the emergence and evolution of life (planetary habitability). 
We study here the role of rotation in the evolution of dynamo-generated magnetic fields in massive Earth- 
like planets, Super Earths (1-10 Mm ) . Usi ng the most recent thermal evolution models of Super Earths 
(|Gaidos et ali I2OIOI: iTachinami et "oD. I2011I) and updated scaling laws for convection-driven dynamos, we 
predict the evolution of the local Rossby number. This quantity is one of the proxies for core magnetic field 
regime, i.e. non-reversing dipolar, reversing dipolar and multipolar. We study the dependence of the local 
Rossby number and hence the core magnetic field regime on planetary mass and rotation rate. Previous 
works have focused only on the evolution of core magnetic fields assuming rapidly rotating planets, i.e. 
planets in the dipolar regime. In this work we go further, including the effects of rotation in the evolution 
of planetary magnetic field regime and obtaining global constraints to the existence of intense protective 
magnetic fields in rapidly and slowly rotating Super Earths. We find that the emergence and continued 
existence of a protective planetary magnetic field is not only a function of planetary mass but also depend 
on rotation rate. Low- mass Super Earths {M < 2 M^) develop intense surface magnetic fields but their 
lifetimes will be limited to 2-4 Gyrs for rotational periods larger than 1-4 days. On the other hand and also 
in the case of slowly rotating planets, more massive Super Earths (M > 2Mq) have weak magnetic fields 
but their dipoles will last longer. Finally we analyze tidally locked Super Earths inside and outside the 
habitable zone of GKM stars. Using the results obtained here we develop a classification of Super Earths 
based on the rotation rate and according to the evolving properties of dynamo-generated planetary magnetic 
fields. 
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1. Introduction 

The number of known exoplanets in the mass range between 1 and 10 Mm is growing 
will be called "Super Earths" or SEs following the classification bv lValencia et all 



hereafter the se objects 
20061 l20n7ahV At the 



time of writing, there ar e almost 46 confirmed planets in this mass range ^ Rivera et a/.l. 12005 :lBeaulieu et al. 
20061: lUdrv et all l2007t iMavor and Udrvl . l2008l iRibas et ali l2008l lOueloz et all 120091: ISonfils et ali I20T1 



Lissauer et al. . 20 111) and more than a few hundred SEs candidates are awaiting further analysis and con- 



firmation (jBorucki et a/.l . 120111 ). These discov eries have increased the inter es t to model and understa i id the 
geophysical properties of this type of planets ( Valencia et al. . 2006 . 2007al lb: Valencia and O'ConnellL 20091 
Seager et all . 2007 : Kalteneggei . 20101 : Korenagall2010 ). The habitability of SEs, in particular those similar 



in composition and structure to the Earth, is an interesting topic in the field and several theoretical works 
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have paid special attention to this particular aspect of SEs properties (lGriel3meier et al. . 2005 . I2OO9I I2OIOI: 



Selsis et al.. 2007: van Thienen et al.. 2007: von Bloh et al. 



20071 iLammer et al\ . miij \. 



Models of th e i nterior structu r e of SEs have been extens i vely developed over the last 5 years ([Valencia et al. 



2006 




2008 





Baraffe et all 120081 : iGrasset et al. 



a,ger et al. . 
%l\. I2OO9I) . 



20061 l2007allbl: iFortnev et al\. 120071: ISeaeer et ali 120071: ISelsis et al\ . 120071: ISotin et aU . 120071: lAdams et al. 



Although there are still open issues to be addressed, 
these models are giving us an understanding of global properties such as the mass-radius relationship 
and its dependence with planetary compositio n, as well as different geophysical phenomena such as man- 
tle convection, degassing and plate tectonics (OlsonL 120071 : IPaouc and Davie^ 2008 : Valencia et al. . 2007 : 
Valencia and O 'ConneliTiooi iKorenagal [2010h . Recently sever al authors have stud i ed in detail the thermal 



evolution and magnetic fi eld properties of this type of planets ( Gaidos et al. . 2010t Tachinami et al. . 2011 



Driscoh and Olsonl . 120111 ) 



Planetary magnetic fields would likely play a role in planetary ha bitability (|von Bloh et am2007l : lGriefimeier et al 



2005ll2009l[2oToHvan Thienen et a/.l . l2007HLammer et al\m\^ . Understanding the conditions for the emer 



gence and long term evolution of a protective planetary magnetic field (hereafter PMF) is crucial to evaluate 
the complex conditions for habitability of SEs. The same conditions could also be applied to evaluate the 
habitability of exomoons around extrasolar giant planets ( Kaltenegger et all . l2010t ). 

The curre nt understanding of PMF emergence and e volution in SEs arises from thermal evolution models for 

the Earth ()Stevensonl . 120031: llabrossi l2003ll2007al lbl: iNimmol . bOOoHA^bert et aLl . boOOHBreuer et aLl.l2010l) 

and scaling laws for convection-driven dynamos obtained from extensive numerical simulations ([Christensen and Auberd . 
20061: lOlson an d Christensenl . l200d lAubert et aU . l2009t IChristensen et all 120091: IChristensenl . l2010l ). Two 



recent works studi ed the problem of PMF evolution in SEs by developing detailed models of planetary 



thermal evolution ( Gaidos et al. . 2010l : Tachinami et al . 2011). Both wor ks have paid special attention to 



different but complementary aspects of the problem. Icaidos et all ( 2010[ ) uses a model of the structure of 
the planetary c ore and its thermal evo lution (hereafter the Core Thermal Evolution or CTE model). On 
the other hand [Tachinami et al. ( 2011 ) uses the Mixing Length Theory adapted to planetary conditions 
to model mantle convection with a detailed treatment of its rheological properties (hereafter the Mantle 
based Thermal Evolution model or MTE model) . Although therm al evolution models for the Earth, other 
terrestrial planets and even SEs have been developed in t he past ([Stevenson , l2003t iLabrossl I2003L fioOTat 
Papuc and Daviesl . [200llNimmol . [2009HBreuer et alihOld ). the CTE and MTE models give the first detailed 



description aimed at studying dynamo-generated magnetic fields of extrasolar terrestrial planets. 
We use the results of the CTE and MTE models to study the role of rotation in the evolution of PMF in SEs. 
We focus on the evolution of the regime of the core magnetic field (CMF) which can be broadly classified 
as non-reversing dipolar, reversing dipolar and multipolar. For this purpose we compute the local Rossby 
number, one of the proxies for CMF regime, as a function of the rotation period and planetary mass. Using 
this property we predict the long term evolution of the surface PMF in rapidly and slowly rotating planets. 
In section [2] we summarize the most important results of the CTE and MTE models. Section |3] presents 
the scaling laws for convection-driven dynamos used to predict the properties of the CMF. In section 2] we 
present a general procedure to compute the CMF intensity in the dipolar and multipolar regimes including 
an implicit dependence on rotation rate. In section[5]we present the results of applying the procedure devised 
here to predict the maximum dipolar component of the field in SEs with different periods of rotation and 
thermal histories as predicted by the CTE and MTE models. Section [6] is devoted to discuss the limitations 
of our procedure and the implications of our results. A summary, concluding remarks and future prospects 
are presented in section [T] For reference a list of the symbols and the physical quantities used in this work 
are presented in Table [T] 
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2. Thermal evolution models of SEs 



2.1. Planetary thermal evolution 

The evolution and long-term survival of a dynamo-g enerated PMfH strongly depends on the thermal history 
of the planet. The energy sources and the amount of energy available for dynamo action change in time as 
the planet cools and, in some cases, its core solidifies. 

Before we present the predicted properties of planetary dynamos and the role of rotation in the evolution 
of those properties, we present the most recent advances in the study of thermal evolution of SEs. Ther- 
mal evolution models are the starting point of any dynamo evolutionary model. The onset of dynamo 
action in the core requires va lues of the magnetic Reynolds number beyond a given threshold Rem > 40 
( Olson and ChristensenL [2006(1 . The strength of the dynamo and hence the intensity of the gen e rated mag- 
netic field are essentially set by the available convective power Qconv I Olson and ChristensenL l2006l) . An 
important role is also played by the size of the convective region that in core dynamos is given by the dif- 
ference between the radius of the outer core and the radius of th e solid inner core. In the thermal evolution 
models studied here ( Gaidos et a/.l . [2010l ; lTachinami et al . 2011 ) Qconv is estimated in two different ways. In 
the CTE model this quantity is obtained by solving first the thermal equilibrium equations for the core and 
mantle and then computing the ohmic dissipation from the balance between entropy sources (radioactive 
decay, secular cooling, sensible heat, latent heat and ot her sources of buoyancy) and sinks (heat conduc- 
tion, ohmic dissipat ion and other sources of dissipation) ( Lister . 2003t Labrosse . 2003 . 2007a : iNimmol . 20091 : 



Gaidos et al. . 2010l ). In the MTE model, Qconv is estimated solving first the thermal transport equations in 



the mantle and the core and computing the difference between the heat that come s out of the core-mantl e 



boundary (CMB) and the amount of energy transported by conduction through it (jTachinami et aLl . l201l[ ). 
Important differences between the CTE and MTE models arise from the way Qconv is estimated (see section 
In the following sections we discuss in detail some of the specific features of the CTE and MTE models. 



2.2. The CTE model 

Gaidos et all ( 2010l ) (CTE model) solved the entropy equilibrium equations in the core assuming adequate 



initial and boundary conditions and a parametrized description of the density and temperature profiles. 
They used a basic model of the heat transport inside the mantle to compute the total planetary energy 
balance and calculate the relevant quantities required to describe the thermal evolution of planets between 
1 and 4.8 M^. Assuming rapidly rotating planets, the y also predicted the evol u tion o f the PMF intensity 
fro m the dynamo scali ng la ws originally develo ped by Christensen and Aubertj ( 20061 ) and later improved 
bv lAubert et al. (l2009l ) and ChristensenI ( 2010l ). In Table [5] we summarize some of the results presented in 
Fig. 8 of iGaidos et a,l\ l{20m . 

In the CTE model the interaction between the core thermal structure, the existence and growth of a solid 
inner core, the mantle properties, the onset of plate tectonics and the surface temperature, determine the 
existence and survival of an intense PMF in SEs. A robust upper bound, close to 2 (low-mass SEs) to 
have strong protective PMF is obtained with this model. Planets with masses larger than this limit will not 
cool enough to develop a solid inner core, a condition that considerably boosts the magnetic field intensity in 
lighter and cooler planets. As a consequence, massive SEs (M > 2M(^) develop a decaying weak magnetic 
field. In this case the dynamo is shut down early when the convective power falls below the condition for 
dynamo action. These results depend strongly on the existence or not of mobile lids (plate tectonics, PT). 



2.3. The MTE model 

iTac hinami et al. ( 201 if ) (MTE model) pays special attention to the role that the mantle have in the extraction 
of heat from the core, taking into account its particular rheological properties. They use the Mixing Length 
Theory modified for solid planets in order to compute the energy transported by convection in the mantle 
and the energy budget at the CMB for planets with masses between 0.1 and 10 M^. They assume a 



^It is important to stress that we are considering in this work magnetic fields generated by dynamo action in a Hquid metaUic 
core. Other fluid sheUs (Hquid, ice or gaseous mantles) may sustain other type of dynamos out of the scope of this work. 
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planet with a mantle-core mass ratio similar to the Earth's. The MTE model includes parameters like the 
rheological properties of the mantle (especially important is the activation volume V* that determines the 
viscosity dependence on pressure and temperature) and the temperature contrast ATcmb in the boundary 
layer between the convecting mantle and the CMB. After solving the mantle and core coupled thermal 
transport equations, they compute, for different model configurations, the heat flux at the CMB Fqmb, 
the heat conducted along the adiabat Fcond at the CMB, and the solid inner core radius Ric as a function 
of time. To predict the PMF intensity and its lifetime they estimate the convective power as Qconv = 
{FcMB — Fcond) X 47ri?^ whcrc Rc is the radius of the outer core. As in the case of the CTE model they 
assume rapidly rotating planets, i.e. planets with a dynamo operating in the dipolar regime. 
Two important new predictions arise from the MTE model: (1) thermal evolution is affected by the strong 
dependence on pressure and temperature of mantle viscosity and (2) the intensity and lifetime of the PMF 
would strongly depend on the initial temperature profile which is characterized by the parameter ISTcmb, 
the temperature contrast between the core and the lower mantle; it is a property mainly determined by the 
early accretion and differentiation history of the planet. Larger values of the initial temperature contrast 
^TcMB, created for example by a violent accrection history, would favor the appearance of an intense and 
long-lived PMF. They found that by doubling ATcmb from 1000 K to 2000 K the dynamo lifetime in 
massive SEs {M > 2 M^) is almost one order of magnitude larger. It should be stressed that the CTE 
model also studied the effects of different thermal profiles by considering the cases of non-habitable surface 
temperatures. This would be the case of highly irradiated planets (not considered in the MTE model). In 
this case the conclusions drawn by both models about the effect of different temperature profiles on the 
PMF properties, were essentially the same. 

2.4- Comparison of the models 

The results of the CTE and MTE models point broadly in the same direction: SEs with masses less than 
approximately 2 Mq seem to have properties better suited for the development of strong and long-lived 
PMFs. The ~ 2 threshold is the most robust prediction of these models. In the CTE model the reason 
that limits the capability of low-mass SEs to sustain a dynamo is the early cooling and stratification of the 
core. On the other hand, the MTE model predicts a strong effect of the pressure-dependent viscosity in the 
thermal evolution of the planet, which limits the capacity of mantle convection to extract heat efficiently 
from the core, especially in the case of massive planets. 

We want to highlight that the evolution of the PMF properties predicted by the CTE and MTE models 
were obtained for planets with short, but not specified, periods of rotation, i.e. planets with dynamos 
operating in the dipolar regime. For rapidly rotating dynamos, the reference dynamo models used by the 
CTE and MTE models (|Christensen and Aubertl . 120061: lOlson and Christensenl . 120061: lAubert et aH 120091: 



Christensen et al . 2009() predict that the magnetic energy density and hence the magnetic field intensity. 



depends weakly on the rotation rate. The still open question regarding the effect that long periods of 
rotation would have in the PMF properties predicted by these models. Particularly interesting is to study 
the effect that rotation would have in the determination of dynamo regime. It would be equally interesting 
to estimate the value for the rotation rate where the rapid rotation approximation could be used. These 
questions are particularly relevant in the case of tidally locked SEs as will be discussed in section 15.31 
For the purpose of this work a limited subset of model configurations included in the CTE and MTE 
models and highlighted in Table [5] have been selected. The final goal is to study the effect of rotation in the 
evolution of long-lived protective PMF in habitable SEfl The chosen configurations must meet two criteria: 
(1) habitable surface temperatures (288 K in the CTE model, 300 K in the MTE model) and (2) long-lived 
and intense PMF. In the CTE case the latter condition requires the existence of plate tectonics. Since we are 
interested in computing upper bounds to surface PMF intensities, given the same set of planetary properties 
(mass, composition, surface temperature, etc.) the case where plate tectonics arises will be the best suited 
to set these bounds. 



solid planet is cons i dered habitable if the teperature on its surface is enough to maintain liquid water on the planet 
surface (see lKasting et oil lll993h and references there in). 
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For planets with 1 the configuration selected in the CTE and MTE models must reproduce the intensity 
of the present Earth's PMF and be consistent with the thermal evolution of the Earth, i.e. they must predict 
correctly the time of inner core nucleation. 

In order to understand the role of rotation in the PMF evolution in SEs we have to take into account the 
dynamo properties that depend on rotation rate. To achieve this we shortly review the relationship between 
the thermal state of the planet and the properties of the dynamo in the following section. 



3. Properties of convection-driven dynamos 



3.1. Scaling laws 

It has been suggested from numerical dynamo experiments that the global properties of convec tion-driven dy- 

namos can be expressed in terms of simple power- law functions of a modified Rayle igh number (jChristensen and Aubert 

son and Christens"enL 120061 : lAubert et al\ . 120091 : IChristensen et a/.l . l2009t ). The scaling relationships 
found by these works involves the convective power Qconv, the core geometry, i.e. the outer core radius and 
the vertical height of the liquid core D = Rc — Ric {Ric is the inner core radius), the rotation rate f2 and 
other properties of the core. 

The p ower-based scaling laws used in this work are expressed using the recent parametrization bv lAubert et al 



(I2OO9II where the properties of the dynamo are scaled in terms of the adimensional convective power density 
P, 

Qc 



Here and V are the average density and total volume of the convecting region respectively. The use of 
p in the scaling laws, instead of the mass anomaly related Ra. yleigh number Rag, is justified from physical 
and numerical grounds (see section 2.2 in Aubert et al. ( 2009f )). 

Two basic adimensional quant i ties have been us ed to characterize the global properties of a dynamo 
( Christensen and Aubert . 2006t Christensen . 2010l) : the Lorentz number Lo and the local Rossby nu mber 



Roi. Lo is the adimensional magnetic field strength and is defined by Eq. (22) in lAubert et al\ (j2009r) as 



Lo 



Br 



(2) 



where Brms — (1/^) / B^dV, is the rms amplitude of the magnetic field inside the convecting shell with 
volume V, fio is the magnetic permeability. Roi measures the ratio of inertial to Coriolis forces and is defined 
as 



Roi 



nL 



(3) 



where Urms and L ^ D/l are the characteristic convective velocity and length scale respectively (7 is the 
mean spherical harmonic degree of the kinetic energy spectrum) . By definition Roi will be large for vigorous 
convection, small characteristic length scale or slowly rotating dynamos. Conversely rapid rotating dynamos, 
large characteristic length scale or a weak convection will produce small values of Roi . 
Numerical dynamo expe r iments covering a wide range of physical properties and boundary conditions 



iNumerical dynamo expe r iments covermg a wide range or pnysical properties and Doundary conditions 
(IChristensen and Aubertl . l2006t lOlson and Christensenl . l200a lAubert et all l2009t IChristensen ali 120091: 
Christensen . 2010lj have found the following scaling relationships for Lo and Roi (Eqs. (22) and (30) in 
Aubert et a/.ll2009l ): 



J- fl/2 1/3 

Lo = CLofoLmP ^ 



(4) 



Ro\ = 



ROI 



(1 + x) 



CRol P 



/5 



(5) 
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where fohm is the fraction of the available convective power converted to magnetic field and lost by ohmic 
dissipation, x = Ric/Rc, E = v/{Q,D'^) is the (viscous) Ekman number {v is the viscous diffusivity) and 
Pr/Pm — X/ K is the ratio of the Prandtl to magnetic Prandtl number (A and k are the magnetic and 
thermal diffusivities respectively). For simpl icity we have approxiinated t he scaling law exponents to the 
ratio of the smallest integers as suggested by Olson and Christensen ( 2006f). W e have also introduced here 



the modified local Rossby number following the suggestion by Aubert et al. ( 2009t ) 
The values of the constants clo and cnoi are obtained by fitting the results of numerical dynamo experiments 
with different boundary conditions. It should be noted th at although the scaling law for Lo in Eq. ^ has 
been generally tested only for dipolar dynamos, recently IChristensenI ((2OIQ,) has found that the magnetic 
energy ^ Lo^ still follows a 2/3-scaling law for multipolar dynamos. In that case the consta nt cr.n fitting 



the p roperties of mviltipolar dynamos is smaller than in the dipolar case by a factor 0.6 (jChristensen 



tne p : 
20101) 



Particulary interesting is that the scaling law for Lo provides a simple way to scale Brms with p irrespective 
of the dynamo regime. Using the definition in Eq. ^ a nd the scaling l aw in Eq. Q we can write down an 
expression for the magnetic field intensity (see Eq. 22 in Aubert et al. ( 2009f l) 



B„ns = CBflili-pclJ^oY'^^Dp^"' (6) 

Given the fact that by definition p cx 17-1/3 (see Eq. (dl), we find that Brms is almost independeni|f| of 
rotation rate both in the dipolar and multipolar regimes. However, the dipolar component of the field will 
depend on rotation rate through the local Rossby number in the case of reversing dipolar and multipolar 
dynamos (see section H]). 

3.2. Dynamo regimes 

Dynamos could be broadly classified according to the power spectrum of the magnetic field at CMB in two 
groups: dipolar dominated dynamos, i.e. dynamos where the dipolar component dominate over other higher 
order components, and multipolar dynamos that have a flatter spectrum or a weak dipolar contribution. 
Quantitatively the dipolarity of core dynamos is commonly measured using the ratio fdip of the mean dipole 
field strength at the CMB, Bdip, to the rms field strength summed up to the harmonic degree 12 at the 
same surface Bcmb, 

fd^p = (7) 



It ha s been assumed that dipolar dominated dynamo have fdip > 0.35 (|Aubert et all 120091 : IChristensenI . 



2OIOI ). Conversely when fdip < 0.35 the dynamo is classified as multipolar. Another useful quantity to 
measure the degree of dipolarity is the relation bdip between the rms strength of the field in the shell volume 
Brms and the dipolar component at the CMB, 

, Brms 

bdip - = — ■ (oj 

^dip 

Large values of bdip are typically a signature of a multipolar dynamo although in numerical experiments 
dipolar dominated dynamos could also have large bdip values. However, the contrary is not true: low values 
of bdip are only found in dipolar dominated dynamos. Typical values of bdip can be found in the lower panel 
of Fig. ((H). 

For the pres ent conditions of the geodynamo fdipm — 0-63 (Bdip(Q = 0.263 mT and Bcmb® — 0.42 mT 
(IOlsonL[2007h ') placing it in the dipolar regime. However, the value of bdip for the geodynamo is largely uncer- 



tain. Estimates of Brms (required to compute bdip) can be computed applying pro per scaling laws. Using for 



example the Elsasser number criterion an estimate of Brms® ~ 4 mT is obtained ([Roberts and Glatzmaier 



*It should be recalled that the exact value of the exponent of p in the scaling law for Lo is close to 1/3 but not exactly equal 
to this value. 
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20001 : lOlsonl l2007l) . Dividing by Bd^,, ^, = 0.263 w e obta in &dip© ~ 15. On the other hand power based 
scahng laws as thos e used in this work ( Christensen . 2010l ) could be used to predict Brms ^1.5 mT (see e.g. 
Aubert et al. 20091) and therefore bdipi^ ~ 5. Both values are also compatible with those found in numerical 



dipolar dominated dynamos (see lower panel of Fig. [T]). 

Finally dynamos could be stable in the long-term or exhibit a reversing behavior (see Amit et al. ( 2010() and 
references there in). In most cases dipolar dynamos do not reverse and r eversing dynamos are mul t ipolar , 
with very little (if any) overlap between dipolar and reversing dvna mos (iKutzner and Christensen , 2002), 
although some special dynamos are dipole-dominated reversing (e.g. OlsonI 2007 ). Dipolarity and reversals 
define the "dynamo regime" . 

Studying numerical dynamos in a wide range of conditions, IChristensen and Auberti (|2006t ) discovered that 
the local Ro ssby number is a proxy of dynamo regime ( see Fig. [H. Their finding has been confirmed by 
other works ( Olson and Christensenl 2006 Aubert et 

MI, I2OO9I: iDriscon and OlsonI , l2009t ) over a wider range 



of dynamo parameters. lAubert et al\ (|2009[ ) found that dipolar dominated magnet i c field s are generated 



by dynamos with values of i?0;* below a threshold of around 0.1. IPriscoll and OlsonI (|2009l ) have confirmed 
this result but using as a dipolarity proxy the modified Rayleigh number Rag (see Figs. 2 and 3 in 
iDriscoll and Olson ( 20091 )). It should be stressed that although Ro* could not be the only controlling factor 
of dipolarity, most of the available evidence points to this q uantity as a good proxy for dy namo regime. 
The value of Roi for the geodynamo is estimated at 0.08 (iOlson and Christensen . 20061) . For = 0-35 
the value of Roi = 0.07. This places our planet close to the boundary between dipolar and multipolar 
dynamos. The reason the geodynamo is so close to that particular boundary is unknown (see section [6] for 
a discussion). 

Dynamo regimes are thus separated in parameter space by complex boundaries broadly limited by approxi- 
mate values of i?0;*. Non-reversing dipolar dynamos has Ro^ < 0.04 (irrespective of the type of convection 
and boundary conditions). The critical value of Roi for the transition from dipolar to multipolar regime 
is in the range 0.04 < i?0;* < 0.1 depending for example of the type of convection. In numerical dynamo 
experiments, where several types of convection are considered, this region is populated by a mixture of 
non-reversing, reversing dynamos and multipolar dynamos. We denote this interval in Roi the "reversing 
region". Finally, dynamos with Roi > 0.1 are multipolar (see upper panel Fig. [1]). 



4. The role of rotation in the PMF properties 



Using the results of the thermal evolution models and the scaling laws p r esente d in section [ 3. 1[ we could try 
to constrain the expected properties of the PMF in SEs. iGaidos et al. (|2010[) and iTachinami et al\ (|201ll ) 
performed this task assuming rapidly rotating planets and hence dipolar dominated dynamos. This work go 
further by including the effect of rotation into determining the dynamo regime. This section will develop a 
general procedure to estimate the PMF properties in planets with dipolar and multipolar CMF, a condition 
required to estimate the magnetic properties of slowly rotating planets. 



4.1. Rotation and CMF regime 

In order to predict the evolution of CMF regime and to constrain the evolving PMF intensity, is necessary 
to find a general expression for the local Rossby number as a function of time, planetary mass and rotation 
rate. Replacing p as defined in Eq. (|T|) and using E — v/{VlD'^), the scaling law for Roi in Eq. ([S]) can be 
written as: 



Ro*i = CRoI X 



-1/2 



^-1/3^1/3^2/31 X (A/«;)l/5 



(9) 



In this expression Qconv and x = Ric/Rc are provided directly by the thermal evolution model. Using x is 
possible to compute D = i?c(l — x) and V = (47r/3) fli^(l — x^). T h e core radius and average density are 
scaled using Rc cx M^-^^^ and cx M^-^^"^ following IValencia et all ( 2006h . Finally, the rotation rate D, is 
expressed in terms of the period of rotation P using = 2tt/P (see section l473l for further comments). The 
general expression for i?o* as a function of time, mass and rotation is finally 
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i?or(t,M,P) = C [p-i/«i?-"/«] X [qVL(1 - x)-"Hl X')-'"] X P'" (10) 

Here we have separated the quantities that exphcitly depend on planetary mass (first bracket) and those 
that depend explicitly on time and implicitly on mass (second bracket). The dependence on rotation 
has been isolated in the P*"/^ term. The quantity C depends on the core viscous, thermal and magnetic 
diffusivities that we will assume are nearly constant in time and almost independent of planetary mass. 
In our models C has been set imposing that a 1 planet at t — 4.54 Gyr and P = 1 day has a value 
Ro^it = 4MGyrs,M = 1M^,P = lday) = 0.07. 

As can be seen in equation (ITU|) the value of i?0;* strongly depends on the rotation period P. Slowly rotating 
dynamos (large P) will have large values of i?0;* and hence will be multipolar. It is also interesting to note 
that dynamos arising from a completely liquid iron core (x = 0) will have a smaller local Rossby number 
and hence will be dipolar dominated for a wider range of periods of rotation. In this case, however, in the 
absence of compositional convection the magnetic field strength could be much weaker. Additionally, for a 
fixed convective power, more massive planets will have smaller values i?0;* and they tend to have dipolar 
dominated dynamos. 

4-2. Scaling the magnetic field intensity 

Now the problem of estimating the CMF intensity in the dipolar and multipolar regimes is tackled starting 
from the thermal evolution model inputs Qconv and D = Rc{l — x) l^see 13.11 for a D definition). Using the 
definition of bdip (Eq- ®) and the scaling law for Brms (Eq- ©) the dipolar field intensity at the CMB 
could be computed as. 
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Here the definition of p (Eq. ([T])) was used and replaced D and V in terms of Rc and x- It is customary 
to assume that in the case of rapid r otatin g dy namos hdiv 1 and hen ce Bdip ^ Brms- That was the 
approximation used by iGaidos et al. ( 2010l ) and iTachinami et al. ( 2011 ) to estimate the field intensities 
reported in their works. However, in many relevant cases (see section 15. 4p this approximation will not be 
valid and a proper prescription to estimate bdip as a function of the therma l evolution model is required . 
We have studied the values of bdip for a set of numerical dynamo results dChristensen and Aubert . 20061 : 
Aubert et al . l2009t IChristensenl 120101 lU.R. Christensen. personal communicationl l201ll ) and found that 
although the value of this quantity varies in a wide range, both in the case of dipolar and multipolar 
dynamos, it is possible to find a lower bound fe^i™ that is a function of fdip (see Fig. [T]), 
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where Cbdip — 2.5 and a ~ 11/10 are the best fitting parameters. Using this equation an upper bound to the 
dipolar field at the CMB could be computed directly from Eq. ([TT|) . replacing bdip as fo^ip" from Eq. (H! 
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The constant c^o adopts different values according to the dynamo regime: in the multipolar case {fdip ^ 
0.35), clo — 1 and in the dipolar case {fdip > 0.35), clo — 0.6. Therefore, clo depends implicitly on fdip. 
In summary, although fdip and bdip do not have unique values for dynamos in a given regime, they are 
constrained from above and below respectively, providing an interesting opportunity to constrain the dipolar 
component of the CMF. Estimating the local Rossby number of a planetary dynamo with a given rotation 
rate it is possible to compute the maximum value of fdip attainable by dynamos at that i?0;* using numerical 
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dynamo results (upper panel in Fig. [T]). With this quantity in hand and using equation (|13p it is possible 
to calculate the maximum value of the dipolar component of the CMF. 

It is possible to write down this procedure in the form of a simple algorithm. Given a planet with mass Mp 
and period of rotation P, the maximum dipolar component of the CMF at time t could be computed using 
the following procedure: 



Using the thermal evolution model find Qconv{t) and 
Compute Ro*i using Eq. (ITU)) . 

Using Rol compute the maximum value of fdip as given by an envelop to numerical dynamo results 
in the fdip — Ro^ space (see dashed line in in the upper panel of Fig. [T|). 

Using fdip compute the maximum dipolar component of the CMF, B^°^ using Eq. Values 
adopt ed for CT,n are y 1 for dipolar dynamos {fdip > 0.35) and ~ 0.6 for multipolar dynamos {fdip < 
0.35) (|Aubert all l2009t [Christensenl . l2010t ). 



We are assuming fohm ~ 1 which is consistent with the goal to obtain an upper bound to the CMF. 
Using the maximum intensity of the dipolar component of the CMF and assuming a low-conductivity mantle 
it is possible to estimate the magnetic field at the planetary surface. In the case of a dipolar dominated 
CMF the surface field intensity scales simply as (Rc/Rp)^- But, if the CMF has a more complex power 
spectrum, the properties of the magnetic field at the surface are harder to estimate. However, since we are 
inter ested in the protect i ve pro perties of the PMF and they mainly depend on the dipolar component of the 
field (|Stadelmann et am2QlQI) , we can still compute the maximum dipolar component of the PMF, 
using the formula: 
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[Appendix A| shows that this result is valid irrespective of dynamo regime. 



4-.3. Evolution of rotation rate 

As was said in section 13.21 the dependence on rotation was expressed in terms of the period P instead of 
the rotation rate Vl = 27r/P. P is better suited to study cases where rotation and orbital periods are related 
(tidal locking) or in cases where large tidal effec ts from planet-stel lar interaction or due to hypothetical 



moons introduce simple long-term variation of P (jVarga et al\ . 119981 ) 



The long-term evolution of rotation periods in terrestrial planets is a complex subject that depends on many 
different effects ranging from dynamical conditions at formation; catastrophic impacts; interior processes 
changing the distribution of matter to tidal interactio ns with the central star; close orbiting bodies or other 
bodies in the planetary system (see van HooltsI (200^ and references therein). 

To model the long-term variations of the rotation period, it is necessary to considered two simple extreme 
scenarios: (1) Constant period of rotation Po\ this will be the case for tidally locked planets and also for 
planets that have preserved their primordial rotation, e.g. Mars. (2) Linearly increasing rotation period. 
This would be the case for planets affected by strong tidal damping from the central star, a close big moon or 
other planetary system bodies. In the latte r scenario, and follo wing the models used to study the long-term 
variation of the Earth's period of rotation (jVarga et al\ . Il998h . we have assumed a simple linear variation. 



P{t)=Po+Po{t~t,). 



(15) 



Here Pq and Pq are respectively the rotation period at to and its rate of variation. For an Earth-like rotation 
we assumed to = 4.54 Gyrs, Pq = 2 4 h and Pp ~ 1.5 h Gyr~^, valu es which are compatible with a primordial 
rotation of p\t = 0) = P„, « 17 h (jVarga et all Il998l: iDenis et aL ,2011.) . 



9 



5. Results 



We have computed the evolution of i?o^ using Eq. ([TU]) and the values of Qconv{t) and x(i) provided by 
the CTE and MTE models in the cases highlighted in Table [5] The results for different planetary masses 
assuming Po=24 h, both in the case of constant and variable period of rotation, are depicted in Fig. [5] 
There are important differences between the results obtained with the CTE and the MTE models. These 
differences arise from the magnitude and evolution of the convective power density in both models (see Fig. 
121) ■ While in the CTE model the available power comes directly from the entropy dissipation inside the core, 
in the MTE model the amount of energy available for convection is simply bounded by the energy extracted 
by the mantle through the CMB. As a consequence, the power density is almost one order of magnitude 
larger and grows faster with planetary mass in the MTE case than in the CTE model. On the other hand, 
the energy flux through the CMB falls more rapidly in the MTE model than the entropy dissipated inside 
the core in the CTE model. This effect produces a net decrease in the power density at later times, at least 
in the case of low-mass SEs. 

The most noticeable change in the evolution of the i?0;* happens during the so-called "p-rebound" just after 
the start of the inner core nucleation at time tic- In the CTE model, p goes through a sudden and strong 
increase at tic due to the combined effect of a convective power increase, which is the result of the release 
of latent heat and light elements, and the reduction in the vertical height D of the convecting shell. The 
p-rebound in the MTE model is milder and comes mainly from the reduction in D. In the MTE model 
the CMB flux Fcmb that determines the estimated value of Qconv, is not sensitive to the new sources of 
entropy dissipation in the recently formed inner core. Although we are comparing both models in equivalent 
conditions, it is clear that the MTE model lacks very important details that reduce the likelihood of the 
conclusions drawn from the application of this model in our case. 

Given these fundamental differences we will present separate analysis of the predicted properties of the CMF 
for the CTE and MTE models in the following paragraphs. 

5.1. Role of rotation in the CTE model 

The evolution of i?o^ for the CTE model is depicted in the upper panel of Fig. [21 As expected, the local 
Rossby number varies in time in a similar way as p does (see Fig. ^ . Evolution of the period of rotation 
has an important effect, particularly when the convective power flux has dropped at late times. 
In the case of low-mass planets, the p-rebound is responsible for the most important features of the i?0;* 
evolution. The sudden increase in p and the related decrease in D produce an even faster increase in Ro^. 
This behavior has two effects: (1) Ro^ reaches a minimum value at tic, e.g. Ro\^^i„ — 0.02 for M = 1 
and Ro\^^i^ ~ 0.01 for M = 2 M^, regardless of variations in the period of rotation, and (2) the inner core 
nucleation, at certain periods of rotation (see below) will mark the transition from a dipolar to a multipolar 
regime. 

In the case of massive planets where the solid inner core does not appear in the first 10 Gyrs, Ro^ steadily 
decreases until the dynamo finally shuts down. In this case the ability of the planet to create a dipolar 
CMF at large periods of rotation (slowly rotating planets) is constrained by the value of the local Rossby 
number close to the time of dynamo shut down. We found in this case that Ro^ ^^^^ ^ lO^'^, which is again 
independent of variations in the period of rotation. 

Taking into account that Roi oc P^/^ (Eq. (UHl)) and using the critical value Rol = 0.1, at which all 
dynamos have became multipolar, we found that planets with a constant period of rotation P < Pmui — 
Po{0-^ / Ro^ j^in)^^"^ will have a chance to develop a dipolar dominated CMF at some point in the dynamo 
lifetime. In the case of planets with M < 2 Mq the limit is Pmui ~ 4 — 9 days (the lowest value corresponding 
to the lightest planet, 1 and for more massive planets, i.e. M > 2M^, Pmui ^ 30— 50 days. These 
limits set upper bounds of what should be considered rapidly rotating planets, i.e. planets able to develop 
strong and long-lived dipolar magnetic fields. To find more precise limits it is necessary to compute the 
dipolar CMF lifetime, a problem that will be addressed below. 

In order to understand the role of rotation in the survival of a dipolar dominated CMF in the case of planets 
with M < 2Mq, the evolution of Ro^ must be studied for several values of the reference period of rotation 
Po in the case of constant and variable rotation rate. The case for a 1 planet with Pq w 1.4 days (33 
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h) is shown in Fig. 2] In this case the CMF dipolarity is guaranteed only during the first few Gyrs when 
no sohd inner core has been formed and the convective power density is still down (low Roi). Shortly after 
the p-rebound the convective forces increase enough to make the CMF multipolar and the period of dipolar 
dominance ends. The time spent by the dynamo in the dipolar dominated regime is called the dipolar 
lifetime, Tdip- More massive SEs M > 2Mq exhibit a different behavior. For periods of rotation P < Pmui 
the CMF is multipolar at the beginning and turns dipolar at a time called the dipolarity switch time, tsw 
In this case T^ip is computed as the difference between the dynamo lifetime and the tsw 
A plot of Tdip for different planetary mass as a function of the reference period of rotation Pq is depicted 
in Fig. [5j As expected T^ip is equal to the dynamo lifetime in the case of rapidly rotating planets (Pq ~ 1 
day). The critical period of rotation Pc up to which this condition is met, increases with planetary mass for 
M < 2 as a consequence of a later formation of the inner core in more massive planets (see Fig. [5]). In 
the case of planets with M > 2 M^, the critical period Pc is almost the same, Pc — 2 — 2.5 days, a behavior 
explained by the similarity between the evolution of Po^* in this type of planets (see upper panel in Fig. [5]). 
A variable period of rotation decreases Pc for planets below 2 Mq but increases in the case of more massive 
planets. The explanation of this behavior could be found in the upper panel of Fig. [2j Low mass planets 
has a lower Pc if Ro^ is larger at late times and this is exactly the effect of a variable period of rotation. 
Conversely more massive planets has a larger Pc if Ro^ is lower at early times as expected again in the case 
of a variable period of rotation. 

It is also interesting to note that in the case of low mass SEs, an intermediate level of dipolar CMF lifetime 
beyond Pc is found. This intermediate level is a by product of the steep increase in Ro^ around the time of 
inner core nucleation tic. For times t < tic the dipolarity condition is ensured in a wide range of rotation 
rates. This effect explains why the dipolar lifetime at this intermediate level is close to tic- With these results 
it is possible to set a more stringent limit to the periods of rotation required to have long-lived {Tdip ^ 3 
Gyrs) dipolar dominated CMFs. We found that regardless of the mass, SEs with P < 2 — 3 days will meet 
this condition. 

An interestin g and nontriv i al pre dictio n arises from the dipola r lifetime dependence on mass and period. In 
the works bv lCaidos et al. ( 201Clf ) and Tachinami et al. ( 2011 ) low- mass planets were identified as the best 



candidates to have intense, long-lived PMF. This result was a consequence of the favorable conditions that 
the early nucleation of a solid inner core and a lower viscosity mantle have on the determination of the PMF 
properties and its lifetime. This result is still valid here, at least for planets with periods of rotation smaller 
than ^1.5 days. However, when rotation periods are larger than this limit the roles are interchanged: 
massive Super Earths would develop a dipolar CMF for a longer time than the lightest planets. The reason 
for this switch is the strong p-rebound when the inner core starts to nucleate in light SEs. Since the thermal 
evolution of massive planets does not exhibit such rebound, the rotation period for M > 2 could be 
increased even more, before the CMF becomes multipolar. 

Taking into account the evolution of Ro* and the information it could provide on the CMF regime, we 
have computed the intensity of the dipolar component of the PMF, B^^dip following the procedure outlined 
in section 14.21 The result is depicted in Fig. El where w e have additionally included the measured value 
of the Earth magnetic field at m'ese nt times (|Maus et all . [2005) and three recent measured paleomagnetic 
intensities at 3.2 and 3.4 Gyrs ago ( Tarduno et al. . 2010l ) (indicated with an Earth symbol and error bars 
respectively). It is clear in this case that the evolution of the Roi has an important impact on the measured 
field at the planetary surface. In the case of an Earth mass planet with P — 1.5 days (lower panel in Fig. 
ini) the dipolar component of the field reaches a maximum intensity almost 300 Myrs after the initiation of 
the inner core nucleation just to decay again in hundreds of megayears as a consequence of the increase in 
intensity of the convective currents vs. the weak Coriolis force. At 3 Gyrs the CMF becomes fully multipolar 
and the intensity measured at the surface is just the CMF dipolar component. We have assumed that the 
transition from the dipolar to the multipolar regime happens in short times when compared with the thermal 
evolution time scales. 

In order to compare the dipolar PMF intensities of planets with different masses and to understand the 
global role of rotation in the determination of the PMF properties, the dipolar field intensity averaged over 
the dynamo lifetime Tdyn was computed 
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Bavg =< Bs^dip >= T^y^^ J Bs^dipdt (16) 

Although not phenomenologicany relevant this quantity is very useful to compare the PMF for evolving 
dynamos with different planetary parameters. On the other hand the minimum and maximum "histori- 
cal" dipolar field intensities for a given planetary mass and rotation rate, observed in a relevant range of 
parameters, are of the same order of this average and therefore Bavg seems a good first order estimate of 
the evolving PMF intensity. In Fig. [7] we plot Bavg vs. Po for planets with different masses. We found a 
similar behavior of Bavg as that found in Tdip (Fig. [5]). In this case however the differences between the 
field intensities are not as noticeable as those found in the case of the dipolar CMF lifetime. 
Field intensity is not the only observable we can try to estimate with this model. Reversal frequencies could 
also be inferred ( DriscoU and OlsonL l2009h . Although we have not computed this property it should be 



mentioned that the predicted monotonous variation of the local Rossby number could not explain completely 
the variation in reversal frequency observed in the Earth's paleomagnetic field. For example the increase in 
Rol after inner core nucleation could explain the transition from a superchrone to reversals but not from 
reversals to superchrone. 

Displaying the properties of an evolving magnetic field (dipolar lifetime and intensity) for different planetary 
masses and periods of rotation is challenging. To simplify the graphical representation of PMF properties 
and their dependence on mass and rotation rate, we introduce here the Mass-Period diagram (hereafter M-P 
diagram). We have used M-P diagrams in Fig. |S]to depict Tdip and Bavg for the CTE model, comparing 
the cases of constant and variable periods of rotation. The conclusions drawn from Figs. [5] and [7] are well 
illustrated in the M-P diagrams (see Fig. Planets with masses between 1 and 3 and periods P < 3 
days, are the best suited to developing a long-lived dipolar field, and planets with M < 2Mq and P < 1.5 
days develop dipolar field intensities of the same order as the Earth's. Qualitatively the effects of a variable 
period of rotation are negligible. 

5.2. Role of rotation in the MTE model 

The evolution of Ro* for the MTE model is plotted in the lower panel of Fig. [51 In all the cases the Ro^ 
has large values at early times and hence the predicted CMF is always multipolar at the beginning. This 
result is a consequence of the large initial CMB flux Fcmb used in this model to estimate the convective 
power at the core. For all the planetary masses and in the case of rapidly rotating planets, dynamos switch 
from a multipolar to a dipolar regime at what we have called the dipolar switch time tsw In several cases, 
the intensity of the p-rebound and/or the variation of the period of rotation reduce the time spent for the 
dynamo in a dipolar dominated state Tdip. In the analogous to Fig. |S]for the CTE model, we have plotted 
tsw and Tdip in M-P diagrams in Fig. [HI In order to have a protective PMF it is expected that the core CMF 
becomes dipolar as early as possible (small tsw ) ■ Indeed a protective PMF i s more important in the early 
phases of planetary and stellar evolution evolution than at late times (see e.g. Lichtenegger et al. 2010f) . We 



also expect that the duration of the dipolarity phase Tdip be also large. The effects of a variable period of 
rotation are noticeable especially in the case of large mass planets. This fact is the result of the particular 
behavior of the evolving Ro^ . 

For the MTE model the limit for rapidly rotating planets, as estimated by the condition to have long-lived 
dipolar fields, is more stringent than in the CTE model. Only planets with periods shorter than 1 day 
develop intense and long lived dipolar CMF in contrast with approximately 2-3 days limit found in the case 
of the CTE model. 

5.3. Rotation based SEs classification 

The resulting effects that a variation in the rotation period has in the dynamo properties of SEs, suggest 
the possibility to classify SE dynamos in three different groups: rapid, slow and very slow rotators (see Fig. 
[SJ. Rapid rotators are planets with strong dipolar dominated dynamos lasting several to tens of Gyrs. 
This type of planets have dynamos with local Rossby numbers below the critical value, i.e. operating in 
the dipolar regime, for the most part of the dynamo lifetime (50% is the criterion used here). Our planet 
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would belong to this group. Slow rotators are SEs with dipolar dynamos lasting for less than 50% of the 
dynamo lifetime. Dynamos of this type spends most time of their lifetime in the reversing region and a 
non negligible time above the critical value for the multipolar regime. Very slovi? rotators correspond to 
planets that do not develop a dipolar dominated field during their entire dynamo lifetime. Dynamos of this 
type have large super critical local Rossby numbers. 

The range of periods of rotation defining the proposed categories, varies with mass and depends on the 
particular thermal evolution model used to predict the PMF properties. In the case of the CTE model 
we found that rapid rotators have P < 1.5 days, irrespective of planetary mass, slow rotators with masses 
M < 2 have periods of rotation in the range 4 < P < 10 days, while more massive slow rotators will 
have 4 < P < 20 days. Planets with P > 10 — 20 days are very slow rotators. We have plotted these limits 
in Fig. O In the case of the MTE model the limit for rapid rotators are more stringent than in the CTE 
model. This fact is a consequence of the high convective power predicted by this model at early times (see 
lower panel in Figs. [2]and[3]). The maximum rotation rate for rapid rotating planets in this case is not larger 
than 1 day. The category of slow rotators is practically nonexistence in this model. Planets with periods of 
rotation larger than 1 — 1.5 days will be very slow rotators if the thermal evolution model has the features 
predicted by the MTE model. 



5.4- Application to already discovered SEs 

The application of these results to study the development of protective magnetic fields in newly discovered 
SEs, depends on the ability to know or estimate their rotation periods. The most of already discovered 
SEs are tidally locked and hence their rotation periods are also known. According to recent claims, planets 
with masses in the range of interest for this work (M < 2 M^) and orbital periods less than 50 days could 
be common around low-mass stars (^How ard et al . 2010|) . Therefore, tidally locked SEs are the best initial 



target for this kind of analysis. Nevertheless, in the case of unlocked SEs yet to be discovered the feasibility 
to measure the pe riod of rotation using p resent and future observational facilities has also been devised (see 
Ford et all (|200lh and lPalle eroZI (|2008^ and references therein). 

In Fig. [TO] we have plotted the habitable zone (HZ) for GKM stars and the ph ysical parameters of confirmed 
SEs plus the subset of Kepler candidat es inside the HZ dBorucki et al . 2011). W e have also included there, 
contours of equilibrium temperatures (IS elsis et a/1 l2007l: iLammer et all . l20l rf). the outer limits to have 
tidally locked planets in 1 and 3 Gyrs (".Kasting et a/. . Il993f ). and the range of periods defining our categories 
of rapid, slow and very slow rotators in the case of massive SEs assuming the CTE thermal evolution model. 
It is needed to point out that the results have a limited range of applicability when compared with the wide 
range of physical properties founded in this set of discovered SEs. We have assumed from the beginning that 
planets have a composition similar to Earth, habitable tem peratures and mob i le lids . Low density planets as 



GJ 1214b probably covered by thick volatile atmospheres (jNettelmann et a/J . l201lh . could have other ways 



to create intrinsic magnetic fields and the limits to develop protective PMF should be different than those 
found here in the case of core dynamos. Nevertheless, rotation will still play an interesting role in the PMF 
properties of these types of planets and efforts to include this effect should not be underestimated. On the 
other hand the composition of many of these SEs is still unknown (there is a la ck of complete inform ation 



on their masses and radii). SEs with high surface temperatures, e.g. Corot 7b ( Valencia et al. . 2010l ). also 
fall outside the range of applicability of these results. Higher surface temperatures reduce the viscosity 
in the mantle favor ing the extraction of energy from the core and increasing the convective power energy. 
ICaidos et all ( 2010l ) found higher values of the surface PMF in the case of w 0(1) Af^ planets and larger 
dynamos lifetime in more massive SEs, when surface temperatures are increased. In those cases the Poj* will 
also be larger and the critical period of rotation to have multipolar dynamos will be smaller. These results 
suggest that hot SEs would lack of protective magnetic fields even with periods of rotation proper of colder 
slow rotators. 



6. Discussion 



This work relies b asical ly on three hypotheses: (1) the thermal evolution models bv lCaidos et al. ( 2010h and 
Tachinami et all ( 201lh provide global robust features of the thermal evolution in SEs, (2) The scaling laws 
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fitted with numerical dynamo experiments can be extrapolated to regions in the parameter space where real 
planetary dynamos lie, (3) the local Rossby number could be used as a proxy for dyna mo regime. 
It is clear that to test hypot hesis (1) it is necessary to address the open questions left bv lGaidos et al. ( 20101 ) 



and Tachinami et al. ()201lh . It is important to solve a complete model including a rigorous treatment of 
convection in the mantle, as was done in the MTE model, but also taking into account a detailed model of 
the structure and entropy balance in the core as done by the CTE model. Despite the limitations, there 
are two robust predictions from these models that may be confirmed by more complete models or even by 
observations: (1) there is a maximum planetary mass, ~ beyond which conditions to develop strong 

and long-lived PMF decline; (2) the formation of a solid inner core is favored in the case of « 0{1)Mq 
planets. These two features are of fundamental importance to our results. Changes in other outputs of 
the thermal evolution model such as the differences between the convective power for planets with different 
mass or the role that other planetary properties will have in the onset of a dynamo, will not change the 
conclusions of this work. 



Hypo thesis (2) is also a matter of concern in present studies of planetary and stellar dynamos (jChristensenl . 
2OIOI ). Altho ugh the application of numeri cal dynamo-based scaling laws to planetary dynamos has had 
some success (jOlson and Christensen . 2006[) . higher resolution in future numerical experiments aimed at 
exploring a wider region of the parameter space is required to confirm this hypothesis. Further advances 
in the understanding of how non dipolar dynamos behave will also be required to support the procedure 
devised in this work to estimate the PMF intensity in that case. 

Hypothesis (3) relies again on infer ences from extensive nu merical dynamo experiments (see Fig. [1]). Since 
the numerical parametric studies of Sreenivasan and Jones| ( 2006) it is well known that dipolarity decreases 
with the increase of the ratio of inertial to Coriolis forces. Christensen and Aubert (I2OO6I) identified a critical 
value for R o* , this critical behavior has been confirmed by more recent studies p erformed by Aubert et 



(|2009l) . iDriscoU and OlsonI ((2009, ) and IChristensenI (|2010l ). lOriscoll and Olso: 



3^ 



in Raq separating the dipo lar and multipolar regiines (se e Fig. 
Ro* critical value found by IChristensen and Aubert (j2006h and Aubert et 



(12009.) found a threshold 
3b in their paper) consistent with the 
([2009'). Less clear are the 

properties of the dynamos lying close to the boundary between the dipolar an d multipolar regime s . The 
best known planetary dynamo, i.e. the Earth's dynamo, is just right there. Driscoll and Olson ( 2009f ) 
discussed the relationship between the particular reversal history of the Earth and the unknown properties 
of the transitional region between dipolar and multipolar regimes in parametric space. In this case however 
an open question remains: why is the Earth's dynamo so close to this boundary? The scaling law for Rol 
found in Eq. ([TU)) could shed light into this "coincidence problem" . It is noted that the particular thermal 
history of our planet does not affect to a large extent the order of magnitude of Ro^. As shown in Fig. [5] 
and [21 one or two orders of magnitude variation of p are not enough to change the order of magnitude of 
Roi . This quantity is more sensitive to the period of rotation and the core radius of the planet. The period 



of rotation of the Earth has been of the same order since the formation of the planet (|Denis et all 120111) 
and it is close to that of Mars. The core radius is mainly determined by the Fe/Si ratio, a quantity that 
is not "fine-tuned" for the Earth since Venus and Mars has a similar value of this ratio. In conclusion the 
present value of the Earth's dynamo local Rossby number is not just coincidentially close to the boundary 
region since, as argued here, reasonable variations in the key properties of the dynamo will place it not too 
far from this region. 

It should be stressed that after the results presented in this work and tak ing into ac coun t the fact that 



proba bly most of the habitable SEs lie in tidally locked regions (see e.g. iBosd 20061 and Forveille et al 
I2OO9I) where large periods of rotation could be common, future efforts to try to understand the emergence 
and evolution of PMFs in these types of planets, must consider the kind of direct and indirect effects that 
rotation have in the PMF properties as those considered here. In other words the assumption of rapidly 
rotating planets, i.e. P < 2 days, is no longer valid if we want to study tidally locked SEs inside the HZ of 
M-dwarfs. On the other hand assuming that tidally locked planets lacks completely of an intense planetary 
magnetic field is also an oversimplification. As has been shown here there are a range of periods of rotation 
where planets could sustain moderate magnetic fields having large periods of rotation 



Thermal evolution models used in this work assumed planetary properties very similar to Earth. iGaidos et 



( 2OIOI ) studied the impact that several modifications to this reference model will have on their results. An 
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interesting case is that of planets with a different relative core size. They found that an increase in core size 
(a larger Fe to Si ratio) essentially has two effects: (1) an earlier nucleation of the solid inner core and (2) an 
increase of the surface magnetic field intensity. The latter effect is mostly due to the smaller attenuation of 
core field and not to a noticeable modification of the convective power density. A different core size changes 
our results in two important ways: (1) the predicted maximum dipolar component of the PMF is increased 
and (2) the dipolar field lifetime, especially in the case of slow rotators, is decreased. Assuming that a 
different core size does not affect noticeably the convective power, the local Rossby number will be slightly 
modified and hence the general results regarding, for example, the intervals of rotation periods for the new 
categories of rotators, will not be substantially altered. 



7. Summary and conclusions 

We studied the role of rotation in the evolution of dynamo-generated magnetic fields in Super Earths. We 
computed the evolution of the local Rossby number and the volumetric magnetic field strength for core 
dynamos in SEs. For this purpose we used the results of two recently published thermal evolution models 
and scaling laws fitted with numerical dynamo experiments. Assuming that the local Rossby number could 
be used as a proxy to dynamo regime, we estimated the maximum dipolar component of the magnetic field 
at the CMB, and from it an upper bound to the dipolar part of the field at the planetary surface. 
We used two properties to characterize the global magnetic properties of SEs: (1) the average of the surface 
dipolar component of the field, Bavg and (2) the total time Tdip spent by the dynamo in the dipolar 
dominated regime ( revers ing and non reversing). Intense magnetic fields with a strong dipolar component 
( Stadelmann et all . hOld ). are best suited to protect planetary environments from external agents (stellar 



wind and cosmic rays). Therefore large values of Bavg, irrespective of the dynamo regime, are consistent with 
planetary habitability. The long-term preservation of water and other volatiles in a planetary atmosphere 
and the development of life, would require long-lived protective PMFs, i.e. large values of Tdip. Intense and 
protective magnetic fields in the early phases of planetary and stellar evolution will be also suited for the 
preservation of an atmosphere or its volatiles. However a planet that achieves to preserve its atmosphere 
during the harsh conditions of the early active phases of stellar evolution but lacks of a protective magnetic 
field soon after this period will leave emerging forms of life to an integrated effect of galactic and stellar 
cosmic rays induced damages. 

We found that the PMF properties depend strongly on planetary mass and rotation period. Rapidly rotating 
SEs P ~ 0(1) day, with mass ~ 0{1) Mq, have the best potential to develop long-lived and intense PMFs. 
More massive planets develop weaker magnetic fields but they have dipolar dominated dynamos in a slightly 
larger range of rotational periods P ~ 1 — 3 days. SEs with rotation periods larger than 3—10 days 
(depending on their mass) will spend the majority of the dynamo lifetime in a multipolar state. 
In order to summarize our results we have introduced a rotation-based classification. Using the CTE thermal 
evolution model, SEs could be rapid rotators P < 1.5 — 4 days, slow rotators 4 < P < 10 — 20 days and 
otherwise, very slow rotators. Planets in the HZ of low mass stars < 0.6 that will be tidally locked in 
less than 1 Gyr, will fall between the slow and very slow rotator types. Unlocked planets could be any of the 
types described before, according to their primordial period of rotation and the effects that could dampen 
it. 

More theoretical and observational efforts should be undertaken to address the problem of direct or indirect 
detection of PMF around low mass planets. The detection and measurement of such planetary magnetic 
fields will help us to constrain thermal evolution and dynamo models. The role of magnetic fields in planetary 
habita bility is another problem that deserves close attention. Recent works have tackled this problem in 
detail ( Griefimeier et al. . 2010f) but their PMF models are too simplistic. Although the effect of rotation 



rates is considered in those models and they have focused on tidally locked planets, their models do not 
include the effects of thermal evolution on the PMF properties and their treatment of the dependence of 
these properties of the rotation rate is also limited. 
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Symbol 


Meaning 


Notes 




Acronyms 


r^iVlr 


Planetary Magnetic Field- 


Surface magnetic field 




CMF 


Core Magnetic Field 


Core surface magnetic field 






V_^Ui C licilllcLi H/VUlUllUn 


vJcLlQUS z-l ilh, iZUlU) 




iVi ± Hi 


iviaiiLie udbeci ±neiniai H/VoiULiun 


±acmiiami cc tit. ^zuiii 






lifxUilcLUlt: Z/Ullt, 

|J 1 o n o T»ir I y* avr i ft c 






itp 


Radius of the core R — M9>0( M / Mn-.)^-'^'^^ 
Average core density, = 1.1 x 10^(M/ Me)^-^'^^ 


Vm VnlpTipin pi nl (OC\C\(\\ 
Kill. V aiciicia Lit. I zuuu ) 

km. Valencia et al. f2006) 




Pc 


kg m-3. Valencia et al. f2006) 






Rotation rate, period of rotation, T = 27r/ri 


rad s~^, days 






Radius of the sohd inner core v — H,:^/H^ 


km 




D 


Vertical height of the liquid core, D — Rc — Ric 


km 




V 


Volume of the dynamo region, V = 4/37r(_Ri? — ) 

Dynamo Properties 


km^ 




Qconv 


Total convective power 






P 


Total convective power density 


Adimensional 




Lo 


Lorentz number, Lo ~< Emag >^^^ 


Adim.. Christensen and Aubert 


(20061 


Ro 


Rossby number, Ro ^< Ekin >^^^ 


Adim.. Christensen and Aubert 


(20061 


Roi 


Local Rossby number, Roi ~< 1^ Ej^i^ ^^/^ 


Adim., Christensen and Aubert (2006) 


fohm 


Fraction of ohmic dissipation 

Magnetic Field Properties 


Adim.. Christensen and Aubert 


(20061 




rms amplitude of the magnetic field inside the convecting shell 






B dip 


Dipolar component intensity of the CMF 






^ s.dip 


Dipolar component of the PMF, Bg^dip = Bdip{Rcl RpY 






fdip 


Dipolar fraction of the CMF, fdip = Bdip/BcMB 


Adim.. Christensen and Aubert 


(20061 


^dip 


Ratio between the rms strength of the field 

and the dipolar component at the CMB, bdip — Brms/ Bdip 


Adim., Christensen et al. (20091 




^ic 


Starting time for the inner core nucleation 


Gyrs 




Tdip 


Dipolar lifetime 


Gyrs 




^sw 


Dipolarity switch time 


Gyrs 






Dynamo lifetime 


Gyrs 





Table 1: Symbols and quantities used in this work. 



Appendix A. Prom the CMF to the PMF intensity 

Assuming a low conductivity mantle and neglecting other field sources, the magnetic field in the region 
outside the conducting core is derived from a pote ntial given by t he solution of the Laplace equation with 
Neumann's boundary conditions at the CMB (Bac kus et aLl . 119961 ). 

V{rA4>) = ~-T.{-) T.<Yr{OA) (A.1) 

where are the spherical harmonics and a[" are the expansion coefficients. Using the harmonic expansion, 
the field regime (dipolar or multipolar) could be described in terms of the power spectrum Wi (r) that expands 
the energy density at radius r, 

tic 



1=1 1 = 1 
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CTE model {Uc,Bs,Tdyn) 



M (Me) 


Tectonics 


288K 


Ts 

1500K 


IMe 


PT 


2.8, 


90, >10 


1. 7,140, >10 




SL 


6.5, 0, >10 


5.9, 0, 


>10 


1.5 Me 


PT 


4.4, 


20, >10 


2.7, 130, >10 




SL 






6.5, 0, >10 


2 Me 


PT 


6.8, 


20, >10 


4.2, 90, >10 


2.5 Me 


PT 


>10, 20, 7 


>10, 30, 6.4 


3.0 Me 


PT 


>10 


20, 6.5 


>10, 30, 10 


4.0 Me 


PT 


>10 


20, 5.2 


>10, 30, 9 






MTE model {U,,Bs,Tdy„ 


) 




M (Me) 


y* 


lOOOK 


2000K 5000K 


lOOOOK 


1.0 Me 


3 m'^ mol^^ 


4, 80, >20 


6.5, 110, >20 


7.5, 130, >20 


7.5, 130, >20 




10 ni'^ mol^^ 


2.7, 80, 10 


2.8, 80, 10 


2.8, 80, 10 


2.8, 80, 10 


2.0 Me 


3 m'^ mol~^ 


0, 90, >20 


7, 120, >20 


8, 140, >20 


8, 140, >20 




10 mol~^ 


0, 0, 0.5 


14, 100, >20 


14, 100, >20 


14, 100, >20 


5.0 Me 


3 m'^ mol"'^ 


0, 0, 1 


7.5, 130, >20 


11, 160, >20 


11, 160, >20 




10 mol^^ 


0, 0, >20 


>20, 0, >20 


>20, 150, >20 


>20, 150, >20 



Table 2: Summary of results for the evolution of PMF in the CTE and MTE models l|Gaidos et al I I2OIOI : iTachinami et alV 
|2011| '|. For every mass and each pair of independent planetary properties (tectonics and surface temperature, Ts in CTE model, 
activation volume, V* and temperature contrast at CMB, ATcjv/s in MTE model), we present the value of three properties 
of the dynamo and the predicted PMF: ti^ (Gyrs) the time for the starting of the inner core nucleation, Bs{to) (fJ-T) surface 
magnetic field at a reference time taken here as the present age of the Earth, 4.54 Gyrs and T^yn (Gyrs) the lifetime of the 
dynamo. All the values are approximated and have been used to characterize the global conditions to have a protective PMF. 
In the CTE model stagnant lid (SL), as opposed to plate tectonics (PT) configurations, are not able to produce a dynamo for 
masses larger than 1.5 and are not included in the Table. We have highlighted the configurations used in this work to 

study the role of rotation in the PMF evolution. 
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Figure 1: Upper panel: observed values of the f ji„ ra t io in numerica l dyna mo e xperiments. CA06, ALP09 and Cll stand 
for da ta obtained from IChristensen and Auber^ ( 20061) , lAubert et al\ d200Sf) and lU.R. Christensen. personal communication! 
J201lfl respectively. The dashed curve is the maximum value of f^ip attainable at a given value of the local Rossby number. 
The regions corresponding to different dynamo regimes have been schematically depicted using shading vertical bands. Given 
the general complexity of the problem either multipolar dynamos could be found in the shaded region of reversing dipolar 
dynamos or non reversing dynamos {fiUp ~ 1) with certain boundary conditions could be found in the reversing band. Lower 
panel: values of fe^jip and f^ip for the same set of numerical dynamo results used in the upper panel. A correlation between the 
minimum value of fe^jip at a given value of f^np is used in this work to estimate the maximum value of the dipolar component 
of the CMF (see text). The position of the geodynamo is indicated with the symbol ®. In the lower panel the uncertainty in 
ferfip© is indicated with the thick black line. 
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Figure 2: Evolution of the local Rossby number computed from selected results of the CTE (upper panel) and MTE (lower 
panel) models. Solid and dashed lines correspond to the cases of constant and variable period of rotation respectively. Earth- 
like values for the period of rotation were assumed, i.e. Po=24 h and to = 4.54 Gyrs. In the case of a variable period of rotation 
(dashed lines) we have used Po = 1.5 h Gyrs~^. Shaded regions enclose values of Ro* corresponding to multipolar (gray upper 
region), dipolar reversing (yellow middle region) and dipolar non-reversing (green lower region) dynamo regimes. 
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Figure 3: Evolution of the convective power density p in tlie CTE (solid lines) and MTE (dashed lines) models for two different 
planetary masses, p is one order of magnitude larger in the MTE model and falls faster than in the CTE model due to 
differences in the estimations of the convective power. 
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Figure 4: Comparison of the evolution of Ro* for a 1 planet with a period of rotation P 1.4 days and P 1 day in 
the CTE model. The Ro^ in the case of variable period of rotation (dashed and dashed-dotted lines) has been computed for 
two different rates Po assuming the same reference period of rotation Po ~ 1.4 days at t = 4.54 Gyrs. The squares are placed 
at the times where the dynamo becomes multipolar (dotted line at Ro* =0.1 marks the transition to the multipolar regime) 
for the constant and variable period of rotation cases respectively. Notice that the dashed line is below the dashed-dotted line 
before t = 4.54 Gyrs and above it after that time. Shaded regions enclose values of Rol* corresponding to different dynamo 
regimes, multipolar (upper region, Ro^ > 0.1), dipolar reversing (middle region, 0.04 < Ro^ < 0.1) and dipolar non-reversing 
(lower region, Ro* < 0.04). 
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Dipolarity lifetime (CTE) 

■ ■ \iSr.Rot. 




K (days) 

Figure 5: Lifetime of the dipolar dominated CMF in the CTE model obtained from the analysis of the Ro* evolution for 
planets with different masses. Constant (solid line) and variable (dashed lines) periods of rotation have been assumed. The 
gray regions are limited by the maximum rotation periods of low mass planets (M < 2 M^) in the three categories introduced 
in this work (see text). 
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Figure 6: Maximum dipolar component of the surface PMF estimated with the procedure devised in this work and using the 
thermal inputs of the CTE model. In the upper panel the case of rapidly rotating planets, Po = 1 day, is presented. In the lower 
panel the dipolar PMF intensities for a 1 Mg^ planet with different periods of rotation are compared. The value of the ijoj* has 
also been plotted in order to illustrate the effect that a transition between dynamo regimes have in the evolution of the PMF. 
The red continuous, dashed and dotted lines in the ijo* subpanels are the limits between the regimes (dipolar non reversing, 
dipolar reversing and multipolar). In the lower panel, the dashed lines in the magnetic field plot are the lifetime average of the 
maximum dipolar component of the PMF, Bavg (Eq. III6I I'). The value of B^ cLip for the pr esent Earth's magnet ic field (Earth 
symbol ©) and recent paleointensities measurements (error bars) has also been included llTarduno et al ] |2010l : Kaap Valley 
(KVP), Barberton Greenstone Belt (BGB), and Nondweni Greenstone Belt (NGB), dacite localities). 
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Figure 7: Surface dipolar field intensity averaged over the dynamo lifetime, Bayg as defined by Eq. 1II6II . for planets with 
constant (solid) and variable (dashed lines) periods of rotation. 
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M-T: i3,_,.. ,(f(T) - CTE 




/^^ (days) 

Figure 8: Mass- Rotation (M-P) diagrams for the dipolar field lifetime, T^ip and the average surface magnetic field, Bavg in the 
CTE model. In each diagram the case for constant (upper half of each panel) and variable (lower half of each panel) periods 
of rotation have been assumed. Contours of equal values of the quantities represented on each diagram are also included. 
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Figure 9: Mass-Rotation (M-P) diagrams for the dipolar field lifetime, T^ip and the dipolarity switch time ts^ in the MTE 
model. In each diagram the case for constant (upper half of each panel) and variable (lower half of each panel) periods of 
rotation have been assumed. Contours of equal values of the quantities represented on each diagram are also included. 
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Figure 10: Stellar mass vs. orbital period diagram for planets around GKM stars {M < I.OSMq). Circles indicate the 
position of 21 of the known SEs and 4 Kepler candidates orbiting stars in the selected mass range. Each circle has a diameter 
proportional to the minimum mass of the planet. Contours of equal equilibrium temperatures Teq as measured in planetary 
surfaces assu ming an Earth-like bond albedo {A = 0.29) and a redistribution factor of 2 as expected for tidally locked planet 
l lSelsis et al\ .\2007) has been also plott ed. Blue and red l ines are the inner and outer limit of the habitable zone as computed 
using the Venus and Mars criteria in llSelsis et all. |2007| ) . Green line is the 1 AU-equivalent distance where the planet will 
receive the same flux as the Earth llKaltene£gej . l2O10[) . The vertical thick dashed line indicates the maximum distance inside 
which planets will be tidally locked in times less than 1 and 3 Gyrs l lKasting et al. 
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In the case of a dipolar dominated CMF, i.e. a magnetic field where the averaged dipolar component is a 
significative fraction of the total averaged magnetic field strength (see section 15^ . the dipolar contribution 
to the power spectrum Wi{Rc) = wi is generally larger than those from higher order harmonics (see for 



examp le the CMF power spectrum of the Earth and other simulated dynamos in Fig. 4 of lDriscoll and Olson 



([20090). Conversely when the field is multipolar the contribution from wi is of the same order or smaller 



than the higher order harmonics contributions. The dependence on r of Wi{r) changes the power spectrum 
of the field outside the core. If the CMF is dipolar dominated, the surface field will be strongly dipolar. 
However, in the case of a multipolar field the surface field regime will depend on the CMB power spectrum. 
Multipolar fields with a flat spectrum wi ~ wi for I > 1 will be dipolar on the planetary surface. A strongly 
multipolar field with a dominated component of order Imax will exhibit a surface fiat spectrum provided 
wimax — wi{Rc/ Rp)'^~^^"^'^'' ■ Assuming a quadrupole dominated CMF and Rc/Rp ~ 0.3 this conditions 
implies that for W2 ^ lOwi the surface field will still be weakly dipolar. 
The magnetic field outside of the core B{r, 6, (/>) = —WV{r, 9, (f) could also be expanded, 

oo I 
1=1 m=0 

The contribution of order / is Bi{r) — X]L=o^r(^) ~ Bi{Rc){Rc/ry^^ . Using this notation the dipolar 
component of the CMF is written as B^ip —< \Bi{Rc)\ > where the average is computed over the CMB 
surface. On the planetary surface the dipolar component of the PMF will be then given by, 

Bs,dip — Bdip (^J^^ (^-4) 
irrespective of the regime of the surface field. 
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